First cosmological constraints on the Superfluid Chaplygin gas model 
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In this work we set observational constraints of the Superfluid Chaplygin gas model, which gives 
a unified description of the dark sector of the Universe as a Bose-Einstein condensate (BEC) that 
behaves as dark energy (DE) while it is in the ground state and as dark matter (DM) when it is in 
the excited state. We first show and perform the various steps leading to a form of the equations 
suitable for the observational tests to be carried out. Then, by using a Markov Chain Monte Carlo 
(MCMC) code, we constrain the model with a sample of cosmology-independent long gamma-ray 
bursts (LGRBs) calibrated using their Type I Fundamental Plane, as well as the Union2.1 set and 
observational Hubble parameter data. In this analysis, using our cosmological constraints, we sketch 
the effective equation of state parameter and deceleration parameter, and we also obtain the redshift 
of the transition from deceleration to acceleration: z± ■ 

PACS numbers: 98.80.-k, 98.80.Jk, 95.36.+X, 95.35.+d 



I. INTRODUCTION 

Understanding present cosmic acceleration is one of the 
major challenges in modern cosmology. This open prob- 
lem has sparked an intense activity directed to investi- 
gate the possible candidates which fuel this acceleration. 
Typically, this accelerated expansion has been attributed 
to dark energy (DE) which makes up nearly 75% of the 
cosmic substratum. Another unknown component of the 
universe is the dark matter (DM) , the missing mass nec- 
essary to held together galaxy clusters and also needed to 
explain the current large scale structure of the universe. 

In the literature there are a considerable number of 
dark energy models; for a selection of DE models see 
[T]. Between those based on scalar field models we find 
are quintessence [2], k-essence [3], phantom 14], and so 
forth. Another possibility to obtain accelerated expan- 
sion is provided by theories with large extra dimensions 
known as braneworlds [5] . Furthermore, there is a stream 
that considers DM and DE as different manifestations 
of the same component and describes the dark sector 
as some kind of fluid whose physical properties depend 
on the scale: it behaves as DM at high densities and 
transforms into DE at lower ones. Most of these Unified 
Dark Matter models invoke the generalized Chaplygin 
gas (GCG), a perfect fluid characterized by a negative 
pressure which is inversely proportional to the energy 
density, as its paradigmatic example [SHE], but in this 
same vein we also find more phenomenological proposals 

Regarding unified dark matter models, recently in 
[TO] the author has proposed a superfluid Chaplygin gas 
(SCG) as a possible unification of both dark sector com- 
ponents: when in the state of a (background) Bose- 
Einstein condensate (BEC) it behaves like DE, while its 
excited states act as DM. The condensate possesses the 
equation of state of the Chaplygin gas, but the evolution 
of the Universe provided by this matter is different from 
that in the two-component model with the Chaplygin gas 
and cold dark matter (CDM) as well as from the GCG 



model used to unify DE and DM. 

In [TU] it is showed that if v, the parameter that gov- 
erns the DM equation of state, is quite large, the model 
will be consistent with cosmological observations, which 
implies that the pressure of normal component is small, 
and thus the superfluid behaves like CDM. From the su- 
perfluid standpoint it means that the second sound speed 
(name for the effective sound speed in the superfluid the- 
ory) is small, a characteristic that should be taken into 
account for any DM equation of state. On the other 
hand, if k, the parameter interrelating the Bose-Eistein 
condensate (DE) and the excited states (DM), can be 
estimated, it can be used to infer the redshift z t of the 
deceleration/acceleration transition because it is directly 
related to the scale factor. 

The good properties of the SCG model proposed in 
|10) . and the subsequent analysis from a state-finder 
point of view developed in [IT] , justifies to carry out a 
observational confrontation using the most recent compi- 
lation of supernovae la (SNe la) Union2.1 [TJ] to explore 
as first approach its reach and faithfulness. Besides SNe 
la, we use the updated sample of Hubble parameter data 
[T3] and the distance modulus of long Gamma-Ray Bursts 
(LGRBs) 14 . 

The work is structured as follows. In Sec. II, we intro- 
duce the main aspects of the SCG model and discuss how 
the equations are to be rearranged for the cosmological 
constraints analysis. Next, in Sec. Ill, we describe the 
observational data samples chosen for our analysis: SNe 
la, Hubble parameter estimates and GRBs. In Sec. IV 
we discuss the constraints and results obtained, and we 
point out our concluding remarks. 



II. THE SCG MODEL WITH FRW SPACETIME 

The SCG model is based on the action 

s = /(-i4 + £ )^ dV (1) 
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where the Lagrangian C is associated to the generalized 
hydrodynamic pressure function depending on a single 
variable if it is considered a pure condensate and on three 
variables if the gas excitations are included. 

In [lOj it is argued that the two-fluid hydrodynamics 
constitutes an efficient scheme to describe excited states 
basically because it does not depend on details of micro- 
scopic structure of the quantum liquid and then capital- 
izes on effective macroscopic quantities. 

The superfluid and normal component are two inde- 
pendent flows, and due that it is necessary to increase 
the number of independent variables in the generalized 
pressure from one to three [TS] . The latter leads to three 
scalar invariants that can be constructed from the pair of 
independent vectors: the superfluid /i Q and thermal 9 a 
momentum covectors allowing to write the general varia- 
tion of the generalized pressure in a fixed background as 
SP = SC = n a 5[i a + s a 86 a with n a and s a the particle 
number and entropy currents respectively. 

The energy-momentum tensor, in terms of the 4- 
velocities U a and V a , is given by 

T aP = pn c V a Vp + W n U a U p - Pg a p. (2) 

Besides, the particle number current is 

n a = n c V a +n n U a . (3) 

The SCG model considers the cosmic substratum as 
matter in the BEC state with energy-momentum tensor 
and particle number current given by Eqs. ^ and ([3| re- 
spectively, and also assumes the following self-dependent 
condensate ansatz 



p (p, P, 7) = Pc(m) + PnO"> 1) ; 



(4) 



besides that the superfluid background obeys the equa- 
tion of state 



A 

Pc = , 

Pc 

with the excited state described by the relations 



(5) 



As usual, the total energy density satisfies a conserva- 
tion law: 



Ptot + 3-(ptot + Ptot) = 0, 
a 



(8) 



which in this case contains implicitly the interaction be- 
tween DE and DM, expressed by the particle number 
conservation, 



^tot + 3-n to t = 0. 
a 



(9) 



In |10j , from the Fricdmann and conservation law equa- 
tions, the following two dimensionless expressions are ob- 
tained: 



, a 2 1 k 
3(1 + i/)-j = - + - 3 



vp 




(10) 



= 0, 



(11) 



where p = p c /vC4, and k — no/vA, the parameter k 
giving an initial normalized total particle number density. 

Clearly, for any given value of the parameter v we 
would have a different solution for the fluid dynamics; 
we note that in the limit we are interested, v — > 00, the 
quasiparticle pressure, Eq. ([6]), is zero, thus behaving as 
dust-like matter (p n = 0), and Eqs. (JlOjTT) can be solved 
analytically. In this case, Eqs. Q - (111 yield to 



k 2 



Pc = 



(a 3 + k ) 2 



1- 



(12) 



while DM is governed by 



Pn 




(13) 



fj,jn n = (1 - c 2 s )W n , p n = c 2 s W n /(l + v), (6) 

where fi is a chemical potential, 7 a scalar associated 
with the relative motion of the components and c s the 
adiabatic speed of sound. 

For a homogeneous and isotropic spatially flat uni- 
verse, the superfluid and normal velocities are equal and 
thus 7=1. In this case, the Einstein equations reduce 
to the Friedmann ones: 



87rGp t , 



a 
a 



AttG 



(3pt. 



Ptotj 



(7) 



where ptot consists of the condensate density, p c , and the 
normal one, p n = W n — p n , that are interpretable as DE 
and DM densities respectively, and ptot = Pc + Pn- 



The parameter k is clearly associated to the ratio of 
normal and condensate density evaluated at present. 



From Eqs. ( 12 ) and ( 13 ) it can be noted that at the be- 
ginning stage (i.e. for small a) the total energy density is 
approximated by ptot oc a -3 corresponding to a universe 
dominated by dust-like matter. In [10] it is pointed that 
although the same behavior is a feature of Chaplygin gas 
[5] and that even though in this model the condensate 
has the same equation of state, such dependence is due 
to the normal component. 



turns out to be 



With Eqs. ( 12 I and ( 13 ), the first Friedmann equation 



ZH 2 = 8ttG 



Kn 



(a 3 + n y 



(14) 
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and after some simple algebra, we can rewrite it as: 



(i 



where, as usual, fl 



= (1 + ^ )_ A / h2 ^ ( ^^ 2 ■ (15) 



H, 



a-' 



fc 2 + (l + K ) 2 

c0 = 87rG/(3i?op Cj o). Henceforth, 
f2 c = O c just for simplicity. It is straightforward to 
derive that, demanding H(z = 0) — Hq, there is a rela- 
tion between Q c and it being: f2 c (l + kq) = 1. Using 
this condition, we can reduce the dimensionality of our 
problem, thus having only two free parameters, £l c and 
k: 



h 2 _ i k 2 n 2 c + [i + (a 3 - i)n c y 



H 2 



i + k 2 n 2 



(16) 



On the other hand, it is known that the equation of 
state parameter is connected directly to the evolution 
of the energy density, and thus to the expansion of the 
universe. In the case of the SCG model, the knowledge 
of k and kq suffice to determine w as follows: 



Ptot 
Ptot 



k 2 + (a 3 + K ) 2 
1 when a — > oo and w goes to 



Note that w goes to 
when a — ¥ 0. 

Finally, we use the deceleration parameter for the SCG 
model given by 



1 



q = 



1 - 



3a 3 (a 3 + k ) 
k 2 + (a 3 + k ) 2 



(18) 



to obtain the transition redshift z t and thus know the 
prediction of the SCG model for the cosmic acceleration 
start. 



III. OBSERVATIONAL TESTS 

In the following, we use the observational data from 
Union2.1 supernovae data set, the update of Union2 re- 
ported in [12], the observational Hubble parameter data 
given in |13j . the distance modulus of long Gamma-ray 
Bursts reported in |14j to constrain the two free param- 
eters of the model (fi c , k). In order to do that, we use a 
Markov Chain Monte Carlo (MCMC) code to maximize 
the likelihood function C[9i) oc exp[— x 2 (9i)/2] where 
6i is the set of model parameters and the expression 
for x 2 (8i) depends on the dataset used. The MCMC 
methods (completely described in [TBHT5] and references 
therein) are well-established techniques for constraining 
parameters from observational data. To test their con- 
vergence, here we follow the method developed and fully 
described in [T§] . 



universe. Recently, the Supernova Cosmology Project 
(SCP) collaboration released the updated Union2.1 com- 
pilation which consists of 580 SNe la [T2J. The Union2.1 
compilation is the largest published and spectroscopically 
confirmed SNe la sample to date. Constraints from the 
SNe la data can be obtained by fitting the distance mod- 
uli p(z). A distance modulus can be calculated as 



H(zj) = 5 log 10 [.D£ (*_,-, +25 
= 5log 10 [d L (zj,9i)] + po, 



(19) 



where po = 42.38 — 5 log 10 h and h is the Hubble constant 
Hq in units of 100 km s~ 1 Mpc _1 , whereas d,L(zj,0i) is 
the Hubble free luminosity distance defined as 



d L {z,6i) = H D L (z,6i) = (l + z) 



dz' 



o E{z',9i) 



, (20) 



in which E(z, 6j) = H(z, 9i)/H and 9i denotes the vector 
model parameters. 

The x 2 function for the SNe la data is 



xl(Po,l 



580 

E 

3=1 



(p(zj;p ,9i) - pobsjzj)) 2 



(21) 



where the cr 2 corresponds to the error on distance mod- 
ulus for each supernova. The parameter fx Q in the Eq. 
( 19 1 is a nuisance parameter since it encodes the Hubble 



parameter and the absolute magnitude M and because 
of that it has to be marginalized over. Here, we follow 
the alternative method suggested in [2D] [2T] which con- 
sists in maximizing the likelihood by minimizing x 2 with 
respect to /iq- Then one can rewrite the x 2 as 



where 



Cl 



C2 



XsnW = c l - 2c 2P0 + C3M0. 



E"° (p(zj-,po = 0,6*,) - pobsjzj)) 1 



(22) 



(23) 



5MI 



(m(^;a*o = 0,0i)- Pobs(zj)) ^ 



3 = 1 



580 ^ 
C3 = S^)- 



(25) 



The minimization over gives po = C2/C3. So the \ 2 
function takes the form 



A. Supernovae 

SNe la are considered standard candles, so they can 
be used to measure directly the expansion rate of the 



XSN( 



Cl 



C3 



(26) 



Since Xsn = Xsn(Po = 0>^i) ( U P t° a constant), we 
can instead minimize Xsn which is independent of fig. 
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B. Hubble parameter observations 

We use the compilation of Hubble parameter measure- 
ments estimated with the differential evolution of pas- 
sively evolving early-type galaxies as cosmic chronome- 
ters, in the redshift range < z < 1.75 recently updated 
in [13] but first reported in [22 . 

The root idea supporting this approach is the measure- 
ment of the differential age evolution of these chronome- 
ters as a function of redshift, which provides a direct esti- 
mate of the Hubble parameter H(z) = — 1/(1+ z)dz/dt ~ 
— 1/(1 + z)Az/ At. The main strength of this approach is 
the confidence on the measurement of a differential quan- 
tity, Az/Ai, which provides many advantages in mini- 
mizing many common issues and systematic effects. In 
addition, compared with other techniques, this approach 
provides a direct measurement of the Hubble parame- 
ter, and not of its integral, in contrast to SNe la or 
angular / angle-averaged B AO . 

Observed values of H(z) can be used to estimate the 
free parameters of the model and also the best-value for 
Hq by minimizing the quantity 



tage to generate a sample that can be directly combined 
with other data to derive cosmological constraints, it has 
also the drawback that is not completely cosmology in- 
dependent because it requires a cosmological model as 
input. 

More recently, in the direction of the first method to 
calibrate GRBs, in 14J it has been estimated the dis- 
tance modulus to long gamma-ray bursts (LGRBs) using 
the Type I Fundamental Plane, a correlation between the 
spectral peak energy E p , the peak luminosity L p , and the 
luminosity time Tl = Ei SO /L p , where Ei SO is isotropic en- 
ergy. Basically the calibration was done in this way: first, 
the Type I Fundamental Plane of LGRBs was calibrated 
using 8 LGRBs with redshift z < 1.4 and SNe la (Union2) 
in the same redshift range by a local regression method, 
to avoid any assumption on a cosmological model; then 
this calibrated Type I Fundamental Plane was used to 
measure the distance modulus to 9 high-redshift LGRBs 
(see [14] for calibration's details). To derive constraints 
on cosmological parameters of the model concerned we 
use these 9 LGRBs calibrated reported in Table 2 of [2] . 

The x 2 function for the GRBs data is defined by 



^mz^-H^yf^ (27) 



where are the measurement variances. In addition, 
we fix H at the value given in [23], H = 73.8 ±2.4. The 
vector of model parameters, 0j, will be Oi — (fi C5 k). 



LGRBs 



E 

0=1 



H(zj,9i) - /iobs(^) 



/'.; 



(28) 



where (J,(zj) — 5 log 10 [dj, (zj , 0j)] + 25. Notice that we 
have used the standard expression of x 2 given through 
the observed distance moduli just to be consistent with 
the way in which the calibration was done. We have fixed 
H as 73.8 ±2.4 from [231. 



The luminosity of GRBs appears to be correlated with 
their temporal and spectral properties. Although these 
correlations are not yet fully understood from first prin- 
ciples, their existence has naturally suggested the use of 
GRBs as distance indicators offering a possible route to 
probe the expansion history of the Universe. Unlike SNe 
la, GRBs can be detected at redshifts beyond 1.7; how- 
ever, the problem is that GRBs appear not to be stan- 
dard candles and to extract cosmological information it 
is necessary calibrate them for each cosmological model 
tested. 

There have been several efforts to calibrate the cor- 
relations between the luminosity and spectral properties 
of GRBs in a cosmology-independent way and among 
the proposals to use SNe la measurements to calibrate 
them externally are [24ll26| . Indeed, recently it has been 
made an update of the sample of GRBs now consisting 
of 109 objects [27] and new calibrations have been devel- 
oped in order to use GRBs in cosmological tests [271 [28] . 
Nevertheless the cosmological constraints from GRBs are 
sensitive to how these are calibrated and if one calibrate 
them using Type la supernovae tighter constraints are 
obtained than calibrating them internally. In [29 GRBs 
have been calibrated internally, without using any exter- 
nal data sets, and although this method has the advan- 



IV. RESULTS AND CONCLUDING REMARKS 

As we have pointed above, to obtain the best fit param- 
eter values of the SCG model, we applied the maximum 
likelihood method. The parameter space of the model is 
(O c , k) where O c represents the fractional energy density 
for DE, and the value of k is directly related to the scale 
factor corresponding to the transition from the decelera- 
tion to the acceleration epoch. Indeed, in [IT] it has been 
found the restriction k < 0.652 for the dust-like normal 
component, in order to be in agreement with the fact 
that the transition should occur at 0.3 < z t < 1 [30] • So, 
from our analysis we would like to know the prediction 
for the value of k as well as of SI c . 

The best fit parameter values with la error and the 
corresponding values of x 2 are summarized in the first 
three columns of Table I and in Fig. 1 the 68% and 95% 
confidence level contours are displayed. As we can see 
from Fig. 1 and from the values for the X 2 e( j> the best- 
fits occur by using SNe la and the combination of all 
observational dataset. In Fig. 2 we present the evolution 
of the effective equation of state parameter w(z) using 
our best fit. In this case and for current time, it turns 
out to be to = — 0.758±R-S51. 
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n c 


k 


Xred 


fc 


Zt 


SNe la 
Hubble 
GRBs 

SNe la + Hubble + GRBs 


0.7741°:°$ 
0.8451°;°?* 

o.869±g : gs 

fl 771+0.046 
U. i 1 1-0.025 


0.2821°;°^ 
0.1741°;^ 

o.i73t8:18 


0.975 
0.913 
0.846 
0.976 


n 9Qi + 079 

U.ZM1_ 142 

0.1831°;^ 
0.15112;^ 

o.297l° : °*I 


U.DO4_ 200 
n 791 +0.130 
u - ' ■ 4x -0.230 

U.y ( ' -0.277 
n 779+O.I6O 
u - ' ' z -0.135 



TABLE I: Summary of the best estimates of model parameters as well as for the transition redshift 
(deceleration/acceleration) z t . The errors are at 68.3% confidence level. See Figure [IJ 
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FIG. 1: Confidence regions in the (f2 c , k) plane for the 
SCG model. The contours correspond to \o-2o 
confidence regions using: SNe la, the green ones; 
Hubble parameter data, the yellow ones; LGRBs, the 
violet ones; the combination of all observational data, 
the red ones. 
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FIG. 2: Effective equation of state parameter w(z) 
using the best fit for k and fco- The central line 
represents the best fit and the shaded contour 
represents the ler confidence level around the best fit. 




FIG. 3: Deceleration parameter q(z) evolving with 
redshift using the best fit for k and kg. The central line 
represents the best fit and the shaded contour 
represents the la confidence level. 



On the other hand, to obtain the transition redshift, 
it is necessary to know the value for acq, see Eq. (18 1. 



To obtain this quantity, we use the normalization condi- 
tion, from which using the respective best fit value for f2 c 
from each observational dataset, we get the correspond- 
ing value for no. The last two columns of Table I shows 
our results for kq as well as for z t , which is a result of 
solving the equation q(z) = 0. Furthermore, in Fig. 3, 
the evolution of q(z) with z can be seen using once again 
our best fit. 

In our results, the best-fit parameters do not have sym- 
metric errors, thus, we cannot use the standard error 
propagation formula and we have to perform a modifica- 
tion in order to account for these non-gaussianities. Here 
we follow the method developed in |32] which we explain 
in the Appendix. 

Our results can be summarized as follows: 

• First, from the constraints using SNe la and the 
combination of all three observational datasets, we 
obtained an abundance of DE, given through f2 c in 
the SCG model, in agreement with the prediction of 
the ACDM model from WMAP7 data alone, fi A = 
0.73 ±0.03 [SI]. 

• Second, when we use the Hubble parameter obser- 
vations and LGRBs, we get restrictions slightly less 
tighter respect to the ones obtained using only SNe 
la; however, the amount predicted for fl c is not 
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so far from the abundance expected. In the case 
of LGRBs, there probably still exist considerable 
systematic uncertainties which are reflected in the 
value of the xfed an d a l so m the respective confi- 
dence region. 



• Third, if we accept the fact that k should be less 
than 0.652 to guarantee a compatible value with 
Zt, as has been suggested in [TT], we can say that 
ours results are in total agreement with this restric- 
tion. In the ACDM model, the expansion of the 
universe switched from deceleration to acceleration 
at zt = 0.752±0.041 [H] and in the SCG model, ac- 
cording to our results obtained from the joint of all 
three observational datasets, the universe exhibits 
the transition from deceleration to acceleration at 
zt — 0.772to^i35- However, the value for z t as a re- 
sult of GRBs alone, is higher than the results from 
SNe la, Hubble parameter and the combination of 
all of them, but still agrees in la confidence. 



Appendix A: Error propagation 



The constraints on the parameters are given in the 
form 6i^ S fg t ' u , where S6i <u and 5di t d are positive quanti- 
ties. 

The estimated error in a quantity depending on them, 
/(#), will be given by an upper limit 



and a lower one 



\ i=i 



(Al) 



A/* = 



^(mm(AjW,-A/a)) 2 

\ 4=1 



(A2) 



where 



A/iu =/(... 6i + A6 iu ,6 ii+1) , ...)- f{6) (A3) 



• Finally, taking together SNe la, Hubble parameter 
and LGRBs, we obtained as best-fit model param- 
eters (fi c = 0.771+°;°^, k = 0-173±g;ig|). Looking 
at the results shown in Table I, we will see that in 
this best-fit, the SNe la result dominates the con- 
straint on f2 c , whereas the constraint on k is vis- 
ibly dominated by the GRBs result. This can be 
understood if we take in mind that f2 c is a param- 
eter evaluated at the present time, a = 1, and that 
k = no/y/X is evaluated at a — 0, thus GRBs rea- 
sonably being a better tool to constrain k than SNe 
la, whereas SNe la are better than GRBs to con- 
strain f2 c . So, in conclusion, this result is affected 
by the redshifts at which the observational data 
set operate. On the other hand, comparing our 
best-fits with those obtained with SNe la, Hubble 
parameter and LGRBs by themselves, we see that 
the errors of model parameters are shrunk when 
LGRBs is used as a complementary cosmic con- 
straint. Indeed, although the number of LGRBs is 
low, in |14) it is pointed out that this sample yields 
to much stronger constraints due to the control of 
systematic errors which have even been removed 
thus enabling consistence with other probes, e.g. 
Cosmic Microwave Background (CMB) and Baryon 
Acoustic Oscillation (BAO). From our results, we 
can also say that LGRB data are a valuable tool 
in probing models where dark matter and dark en- 
ergy are unified if they are considered together with 
other observational datasets. 



Afu = /(... -A0y, %+!),... (A4) 

This error estimation is based on finite differences. If 
the errors are enough small, it can be refined: A9i tU = 
S9i iU and A9ij = 56ij. In that case one can write 




In Gaussian situations, where A9i tU — A0ij — A6i, 
one gets the standard error propagation formula and 
Af u = Af t . 
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